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ABSTRACT 

The  input  power  of  a  radio  frequency  electrodeless 
discharge  operated  at  4  MHZ  with  a  radial  inflow  of  argon 
at  one  atmosphere  of  chamber  pressure  is  measured. 

It  is  found  that  the  discharge  does  not  necessarily 
operate  at  maximum  efficiency;  that  is,  /2~  a/6  ^  2.5. 
The  plasma  radius  to  skin  depth  ratio  for  the  radio  fre- 
quency system  reported  here  was  found  to  be  12.2,  instead 
of  the  value  /2~  a/ 6  =  2.5,  for  which  maximum  operating 
efficiency  is  expected. 
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e  Electric  permittivity 
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CHAPTER  I 
INTRODUCTION 

The  electrodeless  discharge  transfers  by  induction 
radio  frequency  power  to  gases  contained  within  the  volume 
of  a  solenoid  which  is  carrying  a  radio  frequency  current. 
The  axially  directed  magnetic  field  within  the  solenoid 
induces  an  azimuthal  electric  field  in  accordance  with 
Faraday's  Law.   The  azimuthal  eddy  currents  cause  ioniza- 
tion and  ohmically  heat  the  gas  to  high  temperatures. 
This  process  is  of  interest  to  those  charged  with  the 
development  of  high  enthalpy  wind  tunnels  utilized  in  re- 
entry testing. 

One  of  the  earliest  works  that  is  now  recognized 
as  the  fundamental  work  in  the  field  of  electrodeless 
plasma  production  was  reported  by  J.  J.  Thomson  and  is 
described  by  Hollister  [11   in  an  excellent  historical 
review  of  the  literature,   Thomson  developed  the  trans- 
former model  to  describe  the  discharge  and  attributed  the 
observed  phenomena  solely  to  the  electric  field.   How- 
ever, at  this  time  J.  S.  Townsend  claimed  that  the  elec- 
trodeless discharge  was  magnetic  in  origin  and  published 
his  own  theory  of  the  electrodeless  discharge  [1] ,  which 


Numbers  in  brackets  refer  to  similarly  numbered 
references  in  the  bibliography. 
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disagreed  with  the  Thomson  theory.   McKinnon  [1]  resolved 
the  controversy  by  theoretically  demonstrating  that  two 
kinds  of  electrodeless  discharges  exist:   One  of  electric 
origin  and  one  of  magnetic  origin.   Experimental  observa- 
tions by  Babat  [2]  were  interpreted  as  proof  of  the 
presence  of  both  the  Thomson  electric  discharge  and  the 
Townsend  magnetic  discharge  in  the  same  discharge  simul- 
taneously. 

Radio  frequency  (rf)  electrodeless  discharges  are 
manifest  in  many  visible  forms,  and  those  who  have  seen 
these  discharges  can  testify  to  the  various  geometrical 
configurations  which  they  display,  for  example,  diffuse 
glow,  toroidial  discharge,  cylindrical  plasmoid,  arcing, 
and  others.  However,  all  the  aforementioned  discharges 
have  been  divided  into  two  major  categories: 

(1)  The  discharges  in  which  the  elementary  con- 
ductance currents  are  not  closed  but  are  con- 
tinued by  dielectric  currents  are  termed 
"E-discharges. " 

(2)  The  discharges  in  which  the  conductance 
currents  are  in  the  form  of  closed  curves  are 
termed  "H-discharges . " 

No  standard  terminology  has  been  adopted  and  the  two  modes 
of  discharges  have  been  called  respectively:   Low-power 
and  high-power,  glow  and  arc,  low  pressure  and  high  pres- 
sure, cold  and  thermal,  and  capacitive  and  inductive.   In 


this  study,  only  the  H-discharge  (high-power  mode)  will 
be  of  interest. 

The  electrodeless  discharge  described  herein  is 
a  gaseous  discharge  which  occurs  in  a  flowing  argon  gas 
within  the  volume  of  the  work  coil.   This  coil  carries 
an  alternating  current  of  20  to  80  amperes  at  a  frequency 
of  four  megahertz.   Once  the  discharge  is  ignited  in  the 
high-power  mode,  it  manifests  itself  by  a  bright,  luminous, 
white  cylindrical  discharge,  confined  within  the  volume 
of  the  work  coil.   With  little  or  no  gas  flow  and  at 
pressures  less  than  one  atmosphere,  the  luminous  cylinder 
extends  to  the  wall  of  the  confining  quartz  tube  and  no 
dark  interval  can  be  distinguished  between  the  quartz 
tube  wall  and  the  bright  zone. 

In  order  to  illustrate  some  of  the  visible  features 
of  an  rf  discharge,  photographs  of  the  discharge  at  differ- 
ent flow  rates,  pressure  and  power  levels  are  shown  in 
Figure  1. 

Recent  theoretical  energy  balance  studies  of  power 
transfer  have  been  made  by  Keefer  [3] ,  Eckert  [41 ,  Hollis- 
ter  [1],  and  Sprouse  [5]  which  predict  values  for  power 
transfer  that  are  in  agreement  with  similar  predictions 
from  the  Thomson  theory. 

Although  there  has  been  considerable  theoretical 
work  for  determining  the  rf  power  delivered  to  the  plasma, 
very  few  experimental  results  have  been  reported  in  the 


Figure  1.   Photographs  of  the  Discharge. 


5 
literature.   One  recent  measurement  reported  is  that  of 
Dresvin,  Donskoi,  and  Gol'dfarb  [6]  who  used  calorimetric 
and  spectroscopic  techniques.   The  experimental  measure- 
ment of  input  power  is  important  to  the  determination  of 
the  rate  at  which  the  discharge  converts  magnetic  energy 
into  heat. 

It  is  the  purpose  of  this  work  to  measure  the  in- 
put power  to  the  plasma  by  two  methods.   The  first  method 
uses  the  properties  of  the  discharge  geometry  and  the 
applied  magnetic  field  as  predicted  in  the  Thomson  model 
of  coupling.   The  second  method  measures  the  current  and 
voltage  used  to  produce  the  discharge. 


CHAPTER  II 


THEORY 


I.   THE  THOMSON  MODEL  OF  THE  ELECTRODELESS  DISCHARGE 

In  this  section  the  Thomson  model  will  be  used 
to  develop  a  mathematical  expression  relating  the  plasma 
power,  the  input  power  to  the  plasma  and  the  magnetic 
field.   The  Thomson  model  assumes  the  induction  solenoid 
to  be  long  compared  to  its  diameter,  thereby  neglects 
end  effects.   Furthermore,  it  is  assumed  that  the  plasma 
is  uniform  and  homogeneous  and  that  the  variation  of  the 
fields  is  predominantly  a  function  of  the  radial  coor- 
dinate, so  that  the  magnetic  and  electric  fields  are  given 
by 

B  z    i   Real{B  (p)eia)t)  (1) 

and 

E  *  i^  Real^E(t(p)elwt}  (2) 

Maxwell's  equation  for  the  fields  within  the  discharge 
are 

V  •  E  =  0  (3) 

V  •  B  =  0  (4) 
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V    x    E   =    -    icoB  (5) 


VxB=uJ   +    iu>yD  (6) 


The  current  density  J  is  related  to  the  electric  field  E 
through  Ohm ' s  law 

J  =  oE  (7) 

-> 
The  electric  displacement  D  is  related  to  the  electric 

field  by 


D  =  e  E  (8) 

o 


Using  Equations  (7)  and  (8)  ,  Equation  (6)  becomes 


V  x  b  -  m  oE  +  iwu  e   E  (9) 

o        o  o 


The  displacement  current  within  the  plasma  may  be  neglected 
if 


c  >N  we  (10) 

o 


This  condition  is  seen  to  be  satisfied  since  the  electrical 
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conductivity  is  on  the  order  of  10   mhos,  oj  is  of  the  order 

7     -1  -12 

of  10   sec    and  e   ~  10     farads/meter.   Maxwell's  equa- 
tions now  become 


V  x  e  =  -  iu>B  (11) 


V  x .  b  =  u  oE  (12) 
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When  separated  into  components  in  cylindrical  coor- 
dinates (the  geometry  of  the  discharge/  see  Appendix  A), 
we  obtain  for  the  curl  of  E 

/\     o£  ,     E  ,  ^ 

V  x  e  =  i   ( +   —■£-)  =  -  iooi  B  (13) 

z   ~  z  z 

3  p       0 


and  recognizing  that  B   is  the  only  component  of  B  we 


have 


/-.    d  B        /\ 
V  x  b  =  -  i ,  — -  =    -    i ,u  oE.  (14) 

d  p 


From  Equation  (14) ,  we  see  that 

1   dB 
E _i  (15) 


likewise 


dE,       1   92B 

(j)  _    z 

2 

9p      yQa   3p 


(16) 


When  Equations  115)  and  (16)  are  substituted  into  Equa- 
tion (13),  this  yields 

-1  32B     1   1   3B 

+  iu)B   =  0  (17) 


2  z 

\iqo      9p     p  uqo     9p 


or  finally 


3  B  r)B 

p2  1  +  P  — *-  -  (iwuap2)  B   =  0  (18) 

8p^      3p        °       Z 


This  is  recognized  as  a  Bessel  equation  and  is  readily 
solved  by  making  the  following  standard  substitutions. 
Let 


r  =  p/tou  a  e   '  (19) 


Then  Equation  (18)  reduces  to 

92B  3B      9 

r   1  +  r  — -  -  r  B   =  0  (20) 

3r  dr  Z 


which  has  a  solution  of  the  form  [7] . 


B   =  C,I  (r)  +  C„K  (r)  (21) 

z     1  o        2  o 


where  I  (r)  and  K  (r)  are  Modified  Bessel  Functions  of 

zero  order.   The  constants  C,  and  C~  are  evaluated  by 

noting  that  at  r  =  0,  B  (0)  is  finite,  hence  C~  =  0  and 
3  z  2 

that   at 


r   =   a/oou    a  B       (a/mu    a)    =    B 

o  z  o  a 


thus 


Cn    = 


I       (a/coy    a   e      '     ) 
o  o 
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Thus  Equation  (21)  becomes: 


iTr/4 


I       (p/cou    a      e      '     ) 
B    (p)    =   B^    -2 2 (22) 

t       /     / i7T/4, 

I       (a /cog    a      e  ) 

o  o 


If  we  make  the  following  substitutions 


v  =    p/oou    a      =    /2    p/5  (23) 


x  =   a/oou    a     =    /2    a/6  (24) 

o 


where    6   =    /2/jou    o      is   defined   as    the    skin   depth   ratio   and 
the   plasma    radius    (a)    is    to   be    treated   as    a   parameter. 
Then   Equation    (22)    reduces    to 

-r     ,         iff/4 » 

I    (v   e      '     ) 

B    (v)    -    B      — (25) 

Io(xea    '4) 

The  modified  Bessel  Functions  which  have  complex 
arguments  can  be  reduced  to  "Kelvin  Functions"  which  have 
real  arguments  [7]. 


i  '-/4 
I  (x  e   '  )  =  ber   x  +  ibei  x  (26) 

o  o         o 


I  (v   e1T/4)    =    ber   v  +  ibei  v  (27) 

o  o         o 


Thus  from  Equations  (13),  (26)  and  (27),  the  time  varying 
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magnetic  field  is  written: 


(ber    v   +    ibei    v> 

B(v,t)    =    Ba   2 ±-      elwt  (28) 

(ber   x    +■    ibei   x) 


Since  only  the  real  pact  of  Equation  (28)  has  any  physical 
significance,  Equation  (28)  becomes 

?        2   1/2 
(ber^v>  +  bei  v) 

S  (V,t)  =  i   B cos  (art  +  6  ) 

2        2   1/2 
(ber2x  +  be^x)1'2  (2g) 

where 


,  bei  v  ber  x  -  ber  v  bei  x 

,-IO       O  O       O  /->o\ 

>,  =  tan   (30) 

ber  v  ber  x  +■  bei  v  bei  x 
o     o       o     o 


The  equation  for  the  electric  field  is  obtained  from  the 
magnetic  field  by  performing  the  derivative  indicated  in 
Equation  (15) 

2        2   1/2 
_j.         ^0)6     (ber,  v  +  be.i,v) 

E  (V,t)  =  i.  —  B   — — cos  (wt  +  9  )    (31) 

/2         2  2   1/2 

'ber  x  *  bei  x) 
o       o 


where 


,  bei, v  ber  x  -  ber, v  bei  x 

6_  =  tan    — (32) 

ber,  w  ber  x  +  bei,.;  bei  x 
1     o        1      o 


The  expression  for  the  power  distribution  within  the  plasma 
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can  now  be  developed  with  the  use  of  Equation  (31)  .   From 
Poynting's  theorem  we  have  the  power  transferred  to  the 
plasma 

V2  2T     x 

P  =  o        i   i   E2  pdp  d4)  dz  (33) 

)   i 
-1/2      0   0 

for  the  case  considered  here,  Equation  (33)  becomes 


2    4    2        x  2  2 

2ttSLo'jo    6   B         I       (ber,  u   +    bei,v)\;dv 

P   =   -       I (34) 

4  "  2  2 

0       ^ber    x  +   bei    x) 
o  o 


upon  integration  this  reduces  to 


2 

2TT£uia         ~    1     iber.,  x    ?ber    x   -   bei   x)  bei,x    (ber   x  +  bei    x) 

P  =  b2  -    |—i 2 2_  +  1 9. 2_ 

U      a  x  2        2  2        2 

o  'ber  x  t  bei  x)  (ber  x  +  bei  x) 

o  o       o 

(35) 


Equation  (35)  gives  the  power  in  the  plasma  in  terms  of 
the  discharge  geometry,  the  applied  magnetic  field,  and 
the  plasma  radius  to  skin  depth  ratio  x  =  /2~  a/6, 

II.   PLASMA  RADIUS  TO  SKIN  DEPTH  RATIO 

We  want  to  plot  plasma  power  as  a  function  of 
x  =  /2   a/ 6  because  we  are  interested  in  the  variation  of 
plasma  power  with  plasma  radius  and  skin  depth.   This  is 
accomplished  by  making  the  following  substitution 
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2 

2fr£a)a     „ 

P0  =  Ba  (36) 


Thus  Equation  (35)  can  be  written  as 

P         1     rber,x  (ber   x   -   bei   x)         bei,x    (ber  x  .+   bei   x) 
1  o  o         ,  1  o  o 

P      X 


2        2  2        2 

o      l       (ber  x  t  bei  x)  (ber  x  +  bei  x) 

o       o  o       o 


(37) 


This  ratio,  P/P  ,  the  ratio  of  the  power  in  the  plasma  to 
the  input  power  is  plotted  in  Figure  2.   Note  that  the 
maximum  value  of  the  ratio  occurs  at  x  =  /T~   a/6  =  2.5. 
This  can  be  interpreted  as  the  condition  for  an  impedance 
match  between  the  primary  and  secondary  of  a  transformer; 
hence,  the  condition  for  maximum  operational  efficiency 
of  the  discharge  is 

a/cu  a)  =  2.5  (38) 

o  ' 

Although  this  is  the  optimum  point  of  operation,  this  con- 
dition alone  is  not  sufficient  to  determine  the  size  of 
the  discharge  * 

III.   POWER  MEASUREMENT 

As  Eckert  points  out  [4] ,  the  most  important  experi- 
mental parameter  to  report  for  correlation  between  different 

experimenters  is  the  magnetic  field  B  .   As  it  is  shown  in 

a 

Equation  (36),  the  input  power  to  the  plasma  is  determined 
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by  B  .   B   was  measured  using  the  rf  current-transformer 
a    a 

shown  in  Figure  3.   The  magnetic  field  B  for  a  coil  with 

a 

n  turns  per  unit  length  whose  length  is  long  compared  to 
its  radius  is 


B  =  u  nl  (39) 

a    '  o  ' 


The  RMS  current  I   can  be  obtained  using  the  rf  trans- 

Xj 

former  on  the  straight  portion  of  the  lead  from  the  coil 
to  the  oscillator.   The  relationship  between  the  current 
in  a  straight  wire  and  the  magnetic  field  can  be  obtained 
directly  from  Ampere's  Law.   That  is 

I  «  . —  (j,   i  •  ii  (40) 

°  c 

The  geometry  for  the  field  and  line  of  integration  C 
is  shown  in  Figure  3„   Integration  gives 


i   =  !si  i  eiut  (41) 

2*p   * 


From  Faraday's  induction  law,  the  induced  EMF  of  a  loop 
through  which  this  field  passes  is 


.*• 


9B 


Vind  ■  —    B  •  ndS  -    —  •  ndS  (42) 

ina  dt  J-        M   at 


For  the  loop  shown  in  Figure  3 
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Figure   3.      Diagram  of  Radio  Frequency  Current  Transformer 
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11      =      X  A 


and 

iut 


3B        iutu    I    e" 

n    •    —  = (43) 

3t  2np 


Thus  the  induced  voltage  in  a  loop  normal  to  a  straight 

iut  . 
wire  carrying  a  current  I  e     is 

x   iut   £  Yb 
up  I  e      t    t      dp 

v(t)  =  — —  dz  (44) 

2  TT  J         J  p 

0  Ya 


which  yields 

iut 
up  I  e         yb 

v(t)  =  —2 £  ln(— )  (45) 

2  v  y 

J  a 

Since  we  are  interested  in  the  RMS  value  of  the  induced 
voltage  V.  ,,  then 

T 

(Vind)2  =  i   I   v(t)2dt  (46) 

0 

2tt 
where   T   =   — — ,    the   pericd.      Then 
u 

Yb 
MV   ln   ~  h, 

Vind  =   " <47> 

2  IT 


where  I   is  the  RMS  value  of  the  current  in  the  load  coil. 
Solving  Equation  (4T  for  I   and  substituting  into  Equa- 
tion (39)  gives 
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27m 


B   =  V.  ,  (48) 

a  md 

yb 
hi  I   In  (— } 

ya 


When  Equation  (46)  is  substituted  into  Equation  (34),  we 
obtain  the  input  power  into  the  plasma  as  a  function  of 
the  induced  voltage  V.  ,  and  the  associated  discharge 
geometry. 

o  2  2  2 

8t  n  a         9 

P   =  (V.  ,)z  (49) 

o     n  ind 

o 

The  input  power  to  the  plasma  may  also  be  measured 
by  using  the  voltage  developed  across  the  oscillator 
resonate  circuit  V  and  the  total  rf  grid  current  I  ,  as 
illustrated  in  Figure  4.   Since  the  oscillator  circuit  is 
at  resonance,  the  resonant  circuit  and  the  plasma  appear 
as  a  resistance  to  the  rf  generator.   Hence,  under  these 
conditions  the  rf  input  power  into  the  plasma  is 


P-L  =  IT  V  (50) 


It  should  be  noted  that  due  to  the  high  Q  of  the  work  coil, 
essentially  all  of  the  power  P   is  actually  delivered  to 
the  plasma. 
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Figure  4.   Equivalent  Resonance  Circuit. 


CHAPTER  III 
EXPERIMENT 
I.   PLASMA  DEVICE 

The  cylindrical  plasma  was  contained  in  a  4.5  inch 
diameter,  9  inch  long  quartz  tube  which  was  mounted  verti- 
cally as  shown  in  Figure  5.   The  closed  end  was  fitted  to 
facilitate  water  cooling  of  the  cover  plate  and  to  intro- 
duce the  argon  through  four  1/8  inch  copper  tubes  axially 
directed  and  symmetrically  placed  in  the  quartz  tube.   There 
are  20  holes,  .035  inch  diameter,  drilled  in  each  1/8  inch 
copper  tube  and  the  tubes  are  aligned  so  that  the  initial 
flow  direction  is  radially  inward.   The  other  end  of  the 
quartz  chamber  was  connected  to  a  mechanical  vacuum  pump 
which  was  used  to  pump  the  plasma  device  down  to  approxi- 
mately .55  mm  Hg  prior  to  each  initiation  of  the  discharge. 

The  quartz  tube  and  the  end  cover  plate  were  sepa- 
rated from  the  remaining  portion  of  the  plasma  device  by  an 
O-ring  seal.   This  provided  a.  good  vacuum  seal  when  operat- 
ing at  low  pressures.   The  O-ring  seal  also  provided  elec- 
trical insulation  so  that  the  discharge  was  electrically 
floating.   Additional  electrical  isolation  of  the  plasma  was 
obtained  by  placing  a  one  inch  section  of  plexiglass  between 
the  flanged  mounting  brackets  of  the  existing  section  of 
the  plasma  device  and  the  lower  portion  of  the  quartz  tube. 
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Figure  5.   Plasma  Device. 
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II.   OSCILLATOR 

The  oscillator  which  provided  the  rf  power  for  the 
production  of  the  electrodeless  discharge  is  of  the 
Colpitts  design.   Among  the  variety  of  possible  oscillator 
designs,  the  Colpitts  was  chosen  because  of  its  exception- 
ally good  frequency  stability  under  widely  varing  loaded 
conditions.   A  schematic  diagram  of  the  oscillator  is 
shown  in  Figure  6. 

The  heart  of  the  oscillator  is  a  water  cooled  high 
power  triode  which  is  capable  of  delivering  up  to  50  kilo- 
watts of  rf  power.   The  tuning  capacitors  C   and  C   are 

remote  controlled  motor  driven  capacitors.   C   is  used  as 

P 

a  coarse  adjustment  and  C   is  for  fine  adjustments  of  the 

oscillator.   The  load  coil  of  the  oscillator  L   is  wound 

P 

around  the  quartz  tube  in  which  the  plasma  is  generated. 
The  load  coil  was  constructed  of  1/4  inch  copper 
tubing  water  cooled  for  safety  purposes.   However,  the 
water  was  not  turned  on  for  a  few  runs,  and  no  noticeable 
heating  was  observed.   This  shows  that  very  little  power 
is  dissipated  in  the  wcrk  coil.   To  prevent  arcing  between 
the  coil  and  the  metal  cabinet  containing  the  oscillator, 
the  copper  tubing  was  fed  out  of  the  oscillator  through  a 
one  inch  plexiglass  plate  as  shown  in  Figure  7. 


23 


X  i-i 

U  -H 

O    O  Q, 


,-^mKnmnri 


O  O 


"O    0,0 

ai-H  3> 

cS  <h  OT  ^ 

£   -H 

i   o  a>  <n 

CO    <D    ?     I 

woo 

a 


o 

•H 
-P 

e 

CD 
.C 

O 
CO 

o 
+J 


■H 
O 

w 
O 

>i 

u 

CD 

D1 
CD 
S-l 
fa 

O 
-H 

TJ 

rcj 

as 


CD 

(TV 
-H 

fa 


24 


l/4-in.  Copper  Tubing 


Quartz  Tube 


1-in.  Plexiglass 


Oscillator  Cabinet 


Figure  7.   Copper  Tubing  Exit  and  Plasma. 
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III.   OPERATING  CHARACTERISTICS 

To  produce  the  discharge,  the  test  cell  is  pumped 
down  to  approximately  .001  atmospheres.   The  argon  regula- 
tor valve  is  set  for  the  line  pressure  to  be  in  the  neigh- 
borhood of  150  psi.   Before  igniting  the  plasma,  it  is 
best  to  purge  the  system  by  bleeding  some  argon  through 
the  system. 

The  oscillator  resonant  circuit  is  tuned  by  varying 
C   to  obtain  a  clean  sine  wave  on  the  oscilloscope  with  a 
maximum  voltage  output .   A  pleasantly  pink  glow  should 
appear  within  the  quartz  tube  when  the  oscillator  is 
properly  tuned.   This  discharge  is  of  the  low-power  type 
described  earlier.   The  oscillator  plate  voltage  is  in- 
creased from  2  KV  to  approximately  10  KV.   As  the  plate 
voltage  increases,  the  dc  grid  current  will  increase  from 
50  ma  to  as  much  as  800  ma.   During  this  time,  the  discharge 
will  change  from  pink  to  pink-red,  then  to  a  yellow-white 
color.   At  the  same  time,  the  plasma  is  contracting  away 
from  the  end  corner  plates  of  the  quartz  tube  and  continues 
to  contract  until  the  plasma  is  confined  entirely  within 
the  load  coil.   The  discharge  at  this  point  changes  from 
the  low-power  to  the  high-power  discharge .   If  there  are 
any  leaks  in  the  system;  that  is,  if  air  is  mixed  with  the 
argon,  the  system  will  be  very  difficult  to  ignite  in  the 
high-power  mode . 
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Once  the  discharge  is  in  the  high-power  mode,  the 
gas  flow  and  power  are  increased  simultaneously.   Prior 
to  increasing  the  gas  flow,  the  discharge  has  been  attached 
to  the  wall  of  the  quartz  tube.   The  radial  inflow  of  the 
argon  tends  to  "blow"  the  discharge  away  from  the  walls. 
The  final  stable  configuration  of  the  plasma  is  obtained 
by  careful  adjustment  of  the  flow  rate  and  power  input 
which  result  in  a  "hot  plasmoid"  as  illustrated  in  the 
photographs  presented  in  Figure  1,  page  4. 

Typical  plasma  operating  conditions  are  as  follows: 

(1)  DC  plate  voltage 8  KV. 

(2)  DC  current        .  .  .  .  .1.6  amps. 

(3)  DC  grid  current 050  amps. 

(4)  Chamber  pressure  .....  1  atmos. 

-4 

(5)  Flow  rate        .  ....  5  *  10    lbs/sec. 

At  one  atmosphere  of  pressure  and  with  increased  gas 
flow,  the  discharge  becomes  much  more  intense,  and  vacil- 
lates within  the  discharge  without  touching  the  wall. 
The  vacillating  plasma  can  be  stabilized  by  the  proper 
adjustment  of  input  power  and  gas  flow  rate.   In  addition, 
the  stability  of  the  discharge  is  strongly  dependent  on 
the  manner  in  which  the  flow  is  introduced.   For  example, 
the  discharge  is  quickly  "blown  out"  unless  the  increase 
in  gas  flow  rate  is  accompanied  by  a  corresponding  in- 
crease in  input  power. 
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Reed  [8]  reports  that  the  most  stable  operation  is 

obtained  with  a  tangential  flow  of  the  argon,  which  causes 

the  gas  to  spiral  down  the  walls.   This  tends  to  create 

a  low-pressure  region  in  the  center  of  the  tube,  which 

causes  some  of  the  plasma  to  flow  upstream  countercurrent 

to  the  main  flow.   Contrary  to  Reed's  report,  it  is  found 

in  this  experiment  that  tangential  inflow  tended  to  "blow 

out"  or  extinguish  the  discharge  unless  enormous  amounts 

of  input  power  were  dumped  into  the  plasma. 

The  most  stable  plasma  for  any  given  power  or  gas 

flow  rate  was  found  to  exist  with  a  radial  inflow  of 

argon.   With  radial  inflow,  the  plasma  filled  the  volume 

within  the  solenoid  except  for  a  small  region  near  the 

walls  and  vacillation  was  not  observed.   Radial  inflow 

yields  the  optimum  stability  when  input  power,  gas  flow 

rate,  conductive  and  radiative  power  are  considered.   This 

is  in  agreement  with  recent  theoretical  study  by  Keefer  [3] 

who  predicts  improved  stability  with  radial  inflow. 

IV .   INSTRUMENTATION 

The  instrumentation  required  to  monitor  the  elec- 
trodeless  discharge  can  be  divided  into  two  different 
categories: 

(1)   Those  required  to  monitor  gas  dynamic  effects; 
for  example,  instrumentation  for  mass  flow 
rates,  chamber  pressure  and  line  pressures. 
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(2)   The  electrical  instrumentation  used  to  monitor 
the  rf  energy;  for  example,  the  high  voltage 
rf  and  dc  voltmeters,  ammeters,  and  specially 
designed  instruments  for  particular  plasma 
measurements . 
The  gas  dynamic  system  is  illustrated  in  Figure  8  and  a 
diagram  of  the  electrical  instrumentation  layout  is  shown 
in  Figure  9.   The  primary  instruments  for  monitoring  the 
oscillator  consist  of:   M-l,  the  dc  plate  voltage  of  the 
oscillator  tube;  M-2,  the  dc  grid  current;  and  M-3,  the  dc 
cathode  current.   M-l,  2,  and  3  provided  the  total  dc  power 
into  the  system.   A  detailed  account  of  power  calculations 
are  presented  in  the  theory  section  of  this  work. 

One  of  the  most  important  parameters  of  the  plasma 
is  the  current  through  the  load  coil,  since  this  current 
along  with  geometry  of  the  coil  can  be  used  to  correlate 
data  between  different  investigations.   The  load  coil 
current  was  measured  using  the  rf  current-transformer 
measuring  device  T-l  designed  by  the  author  and  shown  in 
Figure  10.   The  rf  transformer  was  fabricated  from  a  one 
inch  slab  of  plexiglass,  and  mounted  on  the  copper  feed- 
in  line  to  L  ,   The  copper  tubing  was  wrapped  with  three 
or  four  turns  of  teflon  tape  to  prevent  arcing  between  the 
transformer  windings  and  the  copper  line.   The  EMF  gener- 
ated in  T-l  was  fed  out  through  a  shielded  coaxial  cable 
to  the  434  Hewlett-Packard  voltmeter. 
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Figure  9.   Electrical  Instrumentation  Layout, 
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Figure  10.   Radio  Frequency  Current  Transformer 
Measuring  Device. 
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A  high  voltage  vacuum  tube  voltmeter,  M-4,  was  de- 
signed by  the  author  to  measure  the  rf  voltage  across  L  . 
Due  to  the  extremely  high  voltage  and  the  rf  frequencies 
involved,  the  meter  was  designed  to  float  electrically  in 
the  circuit.   This  was  achieved  by  mounting  the  tube,  com- 
ponents, and  meter  in  a  light  aluminum  chasis  which  was 
hung  across  L   with  a  pair  of  sturdy  alligator  clips.   A 
diagram  showing  the  meter  mounting  and  the  circuitory  is 
shown  in  Figure  11. 

The  total  rf  line  current  I-  was  measured  with  a. 
Hewlett-Packard  456  rf  current  probe  M-5.   To  insure  proper 
readings,  the  meter  case  was  attached  to  the  cabinet  con- 
taining the  rf  oscillator.   Since  the  total  RMS  value  of 
the  rf  grid  current  exceeded  the  maximum  current  of  the 
probe,  the  oscillator  grid  line  was  divided  into  four  equal 
paths.   The  measured  value  was  then  multiplied  by  four. 
The  output  of  the  rf  current  probe  M-5  was  fed  out  through 
a  shielded  coaxial  cable  mounted  inside  a  copper  tube. 
The  copper  feed-out  line  was  grounded  to  the  rf  shield 
enclosing  the  room  and  the  oscillator  cabinet.   This  in- 
sured that  the  meter  was  at  rf  ground. 

The  oscillator  frequency  was  monitored  using  an  rf 
frequency  meter  and  a  standard  oscilloscope.   A  schematic 
diagram  of  the  electronic  instrumentation  is  shown  in 
Figure  12. 


33 


Copper 
Tubing 


B-3 


10  Meg 


Figure  11.   High  Voltage  Meter  Mounting  and  Circuit  Schematic 


34 


, — <wtnnn 


o 

•H 
-P 

«3 

6 
CD 
& 

o 
w 

c: 
o 

•H 

+j 
fd 
+j 

a; 

g 

5-1 
■P 
M 

H 


<d 
o 

-H 

■P 
O 
Q) 
rH 
W 


tH 

0) 

U 

•H 

En 


35 
As  a  final  note  concerning  the  electrical  instrumen- 
tation, it  should  be  pointed  out  that  each  measurement  was 
obtained  only  after  tediously  eliminating  parasitic  oscilla- 
tions which  seemed  to  exist  everywhere  in  the  vicinity  of 
the  rf  generator.   All  electrical  instrumentation  near  the 
generator  was  eventually  placed  at  rf  ground  by  the  "cut 
and  try"  method. 

IV .   DATA 

The  input  power  to  the  plasma  was  measured  by  two 
methods,  with  the  generator  in  a  stable  condition;  i.e., 
the  plasma  was  not  vacillating  and  was  not  touching  the 
wall  of  the  quartz  tube  and  operating  at  one  atmosphere 
of  pressure.   The  gas  flow  was  held  constant  throughout 
the  runs . 

The  oscillator  dc  cathode  current  and  plate  voltage 
were  measured  to  obtain  the  total  dc  power  into  the  system. 
Then,  while  monitoring  the  dc  instrumentation  the  induced 
voltage  from  the  rf  transformer  T-l  was  recorded.   After 
recording  the  induced  voltage,  the  meter  output  was  checked 
with  T-l  removed  from  the  circuit  to  ensure  that  the  re- 
corded value  was  not  riding  on  an  unwanted  rf  ground  poten- 
tial . 

Next,  the  total  rf  grid  current,  I  ,  was  measured 
by  connecting  the  output  load  of  the  rf  probe  to  the  434 
Hewlett-Packard  voltmeter.   Again,  as  in  the  previous 
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measurement,  the  circuit  was  checked  for  unwanted  rf  ground 

potentials.   After  recording  the  induced  voltage  (B  )  and 

a 

the  total  rf  grid  current,  I  ,  the  dc  power  level  was 
changed  and  the  procedure  repeated.   Three  different  dc 
power  settings  were  used. 

The  oscillator  efficiency  is  determined  as 


P 
— —  x  100  =  efficiency  (51) 


Pdc 


The  efficiency  of  the  oscillator  was  found  to  be 
77  percent.   This  is  a  typical  value  for  a  class  "C" 
operated  high  power  oscillator.   An  oscillator  efficiency 
of  81  percent  is  reported  by  Dresvin,  Donskoi,  and  Gol ^df arb 
[6]  using  calorimetric  and  spectroscopic  techniques  on  an 
argon  electrodeless  discharge  operated  at  one  atmosphere. 

The  results  of  the  power  measurements  and  the 
oscillator  efficiency  calculations  for  the  three  power 
settings  are  presented  in  Tables  I  through  IV.   Table  I 
gives  the  results  as  predicted  using  Thomson's  theory, 
Equation  (4  9) ,  and  Table  II  gives  the  results  obtained 
from  the  resonant  circuit  current  and  voltage  measurements, 
Equation  (50) . 

The  measured  value  of  /2~  a/6  ranged  from  14.4  with 
a  =  2.26  inch  (at  the  wall)  to  12.2  with  a  =  2,00  inch. 
These  values  of  /2~  a/ 6  were  obtained  by  measuring  the  diam- 
eter of  the  high-power  plasma  discharge  from  photographs 
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TABLE  I 


POWER  MEASUREMENTS  FROM  THOMSON'S  THEORY 


Induced 

Voltage 

(volts) 

Load  Coil 
Current 
(amps) 

Magnetic 

Field 
(Weber/m  ) 

Normalized 
Power 
(watts) 

24 

66.0 

3.32  x  10-4 

4,610 

25 

68.5 

3.53  x  10"4 

5,000 

26 

71.5 

3.63  x  lo"4 

5,500 

±  .5  Volts 


±20  0  Watts 
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TABLE  II 
POWER  MEASUREMENTS  FROM  CIRCUIT  VOLTAGE  AND  CURRENT 

Total                   Coil  Radio  Frequency 

Current                 Voltage  Power 

(amps)                 (volts)  (watts) 

1.8  2,500  4,550 

1.9  2,500  4,900 
2.1                                                  2,500  5,250 


±  .04  Amperes 
b±  50  Volts 


C±  340  Watts 
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TABLE  III 


OSCILLATOR  EFFICIENCY 


Cathode3 

Plateb 

c 
DC 

Normalized 

Current 

Voltage 

Power 

Power 

Efficiency 

(amps) 

(volts) 

(watts) 

(watts) 

(percentage) 

1.25 

4,800 

6,000 

4,600 

77 

1.3 

5,000 

6,500 

5,000 

77 

1.4 

5,100 

7,150 

5,500 

77 

a± 

. 02  Amperes 

b± 

150  Volts 

c±  310  Watts 
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TABLE  IV 


COMPARISON  OF  THOMSON  MODEL  WITH  CIRCUIT  CURRENT 
AND  VOLTAGE  MEASUREMENTS 


DC 
Power 
(watts) 


Plasma 

Rad 

io  Frequency 

Input  Power 

Power 

(watts) 

(watts) 

680 

4,550 

735 

4,900 

790 

5,250 

6,000 
6,500 
7,150 
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of  the  runs  and  from  conductivity  values  obtained  from 

Knopp  [9].   A  value  of  plasma  temperature  of  9/000°K,  [10] 
was  used  to  obtain  the  conductivity  value  of  2,000  mhos/m 
from  the  data  presented  [9] .   Using  the  average  conductivity 
of  2,000  mhos/m  [6]  obtained  from  spectroscopic  measure- 
ments on  a  similarly  operated  argon  discharge,  the  value 
of  /2  a/6  ranged  from  17.2  with  (a)  at  the  wall  to  15.3 
with  a  =  2.00  inch.   Thus,  the  plasma  device  is  operating 
somewhat  below  its  maximum  possible  state.   This  is  in 
disagreement  with  the  position  of  Hollister  [1]  who  states 
that  the  plasma  will  always  adjust  itself  to  a  value  of 
/2  a/6   ~  2.5.   For  the  plasma  to  adjust  to  /2~  a/6  =  2.5, 
with  a  =  2.00  inch,  the  value  of  the  electrical  conduc- 
tivity needed  would  be  on  the  order  of  20,000  mhos/m.   The 
measured  values  of  /2  a/6  are  in  agreement  with  Reefer' s 
theory  [3]  and  Sprouse's  results  [10].   Both  state  that 
the  plasma  is  wall  dominated. 


CHAPTER  IV 

DISCUSSION  OF  RESULTS 

As  can  be  seen  from  Table  IV,  page  40,  the  power 
predicted  by  the  Thomson  model  as  determined  by  Equation 
(35)  is  much  less  than  the  rf  power  available  to  the 
plasma  determined  by  Equation  (50)  from  measurements  of 
circuit  current  and  voltage.   The  cause  of  this  differ- 
ence is  not  understood.   The  assumptions  of  the  Thomson 
model  may  not  be  adequate  for  the  complex  problem  of  rf 
plasma  discharges.   For  example,  the  consequences  of 
neglecting  end  effects  and  of  assuming  the  conductivity 
to  be  constant  both  radially  and  longitudinally  are  not 
known.   The  field  calculations  of  the  Thomson  model  in 
Equation  (35)  depend  upon  a  uniform  slab  of  conductive 
material  within  a  uniform  magnetic  field,  where  as  the 
conductivity  of  the  rf  gaseous  discharge  varies  with  both 
the  radial  and  axial  coordinates  [6]  and  the  uniformity 
of  the  magnetic  field  within  the  plasma  is  not  known. 
Furthermore,  the  calculation  of  power  by  Equation  (35)  re- 
quires the  measurement  of  the  induced  voltage  from  which 
B   and  consequently  P   are  determined.   The  measurement 
of  the  induced  voltage  used  to  determine  B   could  have 
introduced  an  unknown  and  perhaps  large  error,  since  only 
a  one  volt  variation  can  cause  a  difference  of  several 

42 


43 
hundred  watts  of  output  power,  see  Table  I,  page  37.  The 
transformer  T-l  used  in  measuring  B  was  not  calibrated 
because  there  were  no  means  available  for  the  calibration 
on  the  present  experiment.  However,  great  care  was  taken 
to  shield  the  cable  thereby  minimizing  any  induced  ground 
potentials. 

The  fault  may  not  all  be  with  the  Thomson  model. 
These  results  open  to  question  the  validity  of  accepting 
determinations  of  rf  power  by  measuring  the  circuit  cur- 
rent and  voltage  as  being  the  power  that  actually  goes 
into  the  plasma.   It  would  appear  that  further  investiga- 
tion is  needed  to  resolve  the  issue. 

The  meter  movements  used  to  measure  the  currents 
and  voltages  in  Tables  I,  II,  and  III,  pages  37,  38,  and 
39,  were  accurate  to  one  percent  of  the  full  scale  read- 
ing.  The  limits  of  accuracy  as  indicated  in  the  table 
footnotes  were  determined  by  the  delta  method. 


CHAPTER  V 


CONCLUSION 


In  conclusion,  the  input  power  to  the  plasma  was 
measured  by  two  methods:   one  based  on  the  Thomson  model, 
and  the  other  on  the  circuit  current  and  voltage  measure- 
ment.  There  was  a  marked  difference  between  the  results 
of  the  two  methods.   The  reason  for  this  difference  is 
not  fully  understood,  although  several  possible  weaknesses 
of  each  method  are  put  forth  and  discussed.   It  is  be- 
lieved that  further  theoretical  and  experimental  study  of 
the  plasma  power  should  be  undertaken. 
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APPENDIX 


APPENDIX  A 


CYLINDRICAL  COORDINATES 


A  diagram  illustrating  cylindrical  coordinates 
used  in  tne  discharge  is  shown  in  Figure  A-l.   The  curl 
of  a  vector  A  in  cylindrical  coordinates,  is 


V  x  a  = 


•\     /s 


1     PI  A   1 

p  K  (j>   z 


3   3   3 


3p  3(j>  3z 


A   pA,  A 
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Figure  A-l.   Cylindrical  Coordinates 
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